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ABSTRACT 

After recent sensitivity upgrades at the Keck Interferometer (KI), systematic 
interferometric 2 /xm studies of the innermost dust in nearby Seyfert nuclei are 
within observational reach. Here, we present the analysis of new interferometric 
data of NGC 4151, discussed in context of the results from recent dust rever- 
beration, spectro-photometric and interferometric campaigns. The complete 
data set gives a complex picture, in particular the measured visibilities from now 
three different nights appear to be rather insensitive to the variation of the nu- 
clear luminosity. KI data alone indicate two scenarios: the i^'-band emission is 
either dominated to ~ 90% by size scales smaller than 30 mpc, which falls short 
of any dust reverberation measurement in NGC 4151 and of theoretical models 
of circum-nuclear dust distributions. Or contrary, and more likely, the i^-band 
continuum emission is dominated by hot dust (> 1300 K) at linear scales of 
about 50 mpc. The linear size estimate varies by a few tens of percent depending 
on the exact morphology observed. Our interferometric, deprojected centro- 
nuclear dust radius estimate of 55 ± 5 mpc is roughly consistent with the earlier 
published expectations from circum-nuclear, dusty radiative transfer models, and 
spectro-photometric modeling. However, our data do not support the notion 
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that the dust emission size scale follows the nuclear variability of NGC 4151 as 
a -Rdust cx: scaling relation. Instead variable nuclear activity, lagging, and 
variable dust response to illumination changes need to be combined to explain 
the observations. 

Subject headings: galaxies: Seyfert — galaxies: individual (NGC 4151) — galax- 
ies: nuclei — techniques: interferometric 



Introduction 



Nearby Active Galactic Nuclei (AGN) are rosetta stones to understand the astrophysics 
close to an actively accreting super-massive black hole. Sub-parsec resolution is required to 
identify answers to the key questions of AGN physics: (i) How does galactic material flow 
down to the accretion disk to feed luminosities close to the Eddington limit? (ii) Which role 
do play outflows and jets in the energetic interconnection between the AGN and its host? 
and (iii) Are AGN of different luminosity intrinsically similar? Studying the spectral energy 
distribution of quasars and type 1 AGN from the op tical to the far-infrared reveals remark 



ably similar SEP shapes, suggesting similar ph ysics ( lEdelson fc Malkanlll986l : ISanders et al. 



19891 : iKobayashi et al.lll993l : iRiffel et al.ll2006l ). Therefore, resolving astrophysical phenom- 
ena in nearby AGN enables to better understand farther, unresolved nuclei, and to gauge 
correlations between nuclear luminosity and properties of the surrounding circum-nuclear 
environment. Being one of the brightest AGN on the sky, NGC 4151 reveals this unique 
potential of nuclei closer than a few tens of Mpc's in numerous publications. 

In this article, we concentrate on the origin of the near-infrared (NIR) emission of NGC 
4151, and its relation to the nuclear luminosity. NIR-AGN emission is characterized by a 
steep decline from optical wavelengths down to about 1 /zm, and on the l ong wavelength side 
by the quick rise of a separat e IR emission bump peaking at about 3 /zm ( lEdelson fc Malkan 
19861 : IKobayashi et al.l 119931 ) . The origin of this NIR exce ss has been discussed through- 



out the past fou r decades (probably one of the first were iPacholczyk fc WeymannI Il968 



Rees et al.lll969[ ). While it is not in question that the original power originates in the accre- 
tion disk, it is uncertain if the NIR emission excess derives directly from nuclear emission 
processes, or if mostly dust outside the broad line region (BLR) re-radiate s the nuclear 
emission (lEdelson fc MalkanI Il986l : lEdelson et al.l Il988l : iBarvainid Il987l Il992l ) . The proper 
measurement of the nuclear near-infrared emission (and its discrimination from dust-re- 
processed light) is crucial to derive the intrinsic SED of accreting supermassive black holes 
(SMBH) in the centers of galaxies and to understand the radiative processes involved. 
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Even at the relatively detailed linear scale of NGC 4151 (~ 82 mpcmas"^, footnote ^""^ 
in Table |3]), the accretion disk itself and the surrounding BLR are too small to be spa- 
tially resolved by t he current generation of optical or infrared telescopes or interferometers 
( iBentz et al.ll2006l ). However, the circum-nuclear dust distribution is now wit hin reach of 



observations with infra.red telescope arrays a t sub-50 mas angular resolution (ISwain et al. 



20031 : IWittkowski et al.ll2004j : Ijaffe et al.ll2004i ) . Mid- infrared data of th e VLT interferometer 



typically find in Seyfert nuclei 10 um em ission sizes of a few parsecs (ITristram et al.l 12009 



Burtscher et al.ll2009l : iRaban et al.ll2009l ). Dust plays an exceptional role among the circum- 
nuclear components. It is not only assumed to significantly contribute to the near-to-mid- 
infrared continuum radiation of an AGN. Also, a non-spherically symmetric dust distribution 
is widely assumed to explain the type 1/2 dichotomy of AGN emission line spectra. Radia- 
tive transfer models of clumpy distributions of dust clouds extending out of the equatorial 
plane are c urrently the fayored explanation of the steadily increasing amount of observa- 



tional da ta (BarvainisI 



1987 



1992 



Schartmann et al.ll2005t iHonig et al.ll2006l : iNenkova et al. 



2008al jbl: iHonig fc Kishimotol " I2OO9I ). In the following we refer to the dust distribution as 



torus, keeping in mind that its detailed morphology is rather uncertain, might not resemble 
closely a smooth torus. 

It is currently an open debate, how the findings and models at MIR wavelengths connect 
to the inne r most hot dust, which is expected to contribute to the continuum emission at 2 /xm. 
Mor et al.l ( 120091 ) argue for a spatially and chemically distinct inner component, based on 
2 — 35 yum spectra. Also, the interferometric MIR data of NGC 4151 appear to reject 
that the nuclear flux at 10 fim is domi nated by the outer, co oler part of the same structure, 
which dominates the 2 /im continuum (IBurtscher et al.ll2009l ). Such results discourage simple 
attempts to unify the dust emission structures around AGN, and suggest multi-component 
models, or at least a significant change o f dust composition and grain size distribution wi th 
the radial distance to the central engine (jSchartmann et al.ll2005t iKishimoto et al.ll2009al ). 



We report in this article on first results of a new campaign using the Keck Interferome- 
ter (KI) to add observational constraints to the origin of the NIR continuum of NGC 4151. 
The 85 m baseline and the sensitivity of the KI is adequate to resolve the linear distances of 
order 30-200 mpc around NGC 4151 which matches the theoretically calc ulated dust subli- 



mation radii for this sourc e. Recent sensitivity improvements of the KI ( iWizinowich et al. 



20061 : iRadand et al.l [20081 enabled us to repeat the early Swain et al. KI measurement, 
although the variable nucleus of NGC 4151 has been at a significantly fainter state dur- 
ing the time of observation. After a description of the observations (Section [2]), we discuss 



^Current KI performance numbers are given at the website of NASA Exoplanet Science Institute 
(Nexsci): |http : //nexsci . caltech. edu/sof tware/KISupport7v2/v2sensitivity .html) 
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our fi ndings and the implications on the interpr etation of the combined visibihty data set 
(from ISwain et al.l l2003l : iKishimoto et aLl l2009bl . and the here presented observations) in 
Section [31 We concentrate our discussion on issues left open by the previous high resolution 
NIR observations of the AGN: (i) is t he i^-band emissi on rather dominated by unresolved 
accretion disk emission (as favore d by ISwain et al.ll2003l ). or by resolved circum- nuclear dust 
emission (IKishimoto et al.ll2009bl )? (ii) How closely related are ac cretion disk luminosity and 
circum-nuclear dust location in NGC 4151 ( iKoshida et al.l 120091 )? (iii) Is the currently de- 
scribed scatter between direct interferometric and indirect reverberation measurements of 
the torus size a systematic offset reflecting that the respective method probes slightly dif- 
ferent dust reservoirs, or does the scatter rather refl ect the currently achievable accuracy 



of eith er method in deriving the actual dust location ( IKoshida et al.ll2009l : IKishimoto et al. 
2009bl )? The concluding remarks in Section H] summarize our current results and give a 



brief outlook on the project. 



2. Observations 

The 85 m baseline of the Keck Interferometer (KI) is oriented 38° east of north (Colavita 
et al. 2004; Wizinowich et al. 2004). We used the /f'-band (2-2.4 yum) V'^ continuum mode. 
All data shown here are from the white-light channel of the beam combiner, observed at an 
effective wavelength of 2.18 fim. The observations were conducted on Dec. 15, 2008 (UT), 
details appear in T able HI We followed standard observing and data reduction procedures 



(see Section 3.2 in iColavita et al.l l2003l and references therein). KI data are provided to 
the observer in a semi-raw state, and still require estimation and correction for the system 
visibility (or visibility transfer function), which is estimated by observing unresolved cali- 
brators, close in space and time (typical numbers for dose are < 10° and < 0.5 hr). We 
chose calibrators of a i^'-band magnitude similar to the target, to avoid a flux bias in the 
calibrated visibilities {V"^, Table |2]). Our typical scanlength, leading to one visibility point, 
is 200 sec. 

A reliable data calibration is indicated by stable system visibility and the flux ratios 
between both telescope beams over the night. Our observation was conducted at a visible 
seeing of 0.5". To estimate and apply the system visibility, we used the wb/nbCalib-software 
suite by Nexsci. We use the standard deviation over the about 25 individual visibility 
estimations, which make up the 200 s scan, as a flrst estimate of the uncertainties of each data 
point. The resulting statistical noise in the calibrated visibility measurement is 0.01-0.015, 
proving the good observing conditions of the night. This gives the differential intranight 
visibility precision, and can be used to check for visibility slopes with respect to changes 
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Table 1: Observing log. 



Target date (UT) 


R.A." 




proj . B ^ 


calibrators 






[hr] 


[m] 


[m, deg(EofN)] 


[from Tab |2] 


[calib.] ^ 


NGC 4151 Dec.15,2008 


-1.3 


(56.4, 51.6) 


(76.4, 47.5) 


1,2 


0.85 


Dec.15,2008 


-1.1 


(56.0, 54.0) 


(77.8, 46.0) 


1,2 


0.84 



"Hour angle 

^The u, i^coordinates, given here in meters, are the basehne length (B) projected onto the line of sight, u 
points East, v points North. They are equivalent to polar values of the projected baseline, given in the next 
column in meters, and degrees East of North. 

'^The absolute calibration accuracy is at least 0.03, which is shown in Fig. [1] and should be used if comparing 
these values with results from different nights and observing campaigns. The differential intranight visibility 
precision is 0.01 - 0.015, based on the statistical scatter of the measurements. 



Table 2: Properties of the interferometric calibrator stars used for calibration of the instru- 
mental transfer function during data reduction. 



# 


Calibrator 


V/H/K 


Spec. Type '•"^ Ang 


. diameter (mas) ^''^ 


1 


HIP58819 


10.9/8.5/8.4 


KOIII 


0.11 ± 0.03 


2 


HD109691 


8.9/8.9/8.9 


AOV 


0.04 ± 0.02 



"from SiMBAD 

''bolomctric diameter fit from the Nexsci getCal tool. 

Note. — Calibrator stellar diameters significantly smaller than 0.5 mas, are unresolved by the KI. The 
statistical errors given here for the bolometric diameter fit to a black body likely slightly underestimate 
systematic errors of the NIR diameter of stars, but even 0.2 mas uncertainties in the diameter would not 
change the visibility calibration. 
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of the baseline projection. Long-term monitoring of KI data against known binary orbits 
shows that the absolute calibration accuracy of the data from night to night is at least 0.03 
B The difference between the intranight precision, and internight accuracy is visualized in 
Fig. [H We use the 0.03 uncertainty to derive absolute emission sizes (Table [3]), and to check 
for flux-induced size variations (Fig. |2]). 

The NIR flux of NGC 4151 is variable, and a contemporaneous total flux measurement 
is needed in our discussion of the data, to properly compare our visibilities with data from 
other observing epochs. A direct flux cahbration of the total photon counts detected by the 
interferometric camera (called Fatcat) needs to reflect, that Fatcat is fed by single- mode 
fibres. The fibres realize a spatial filtering of ~ 50 mas on sky, to reduce the atmosperic phase 
noise. Since this quasi- aperture is on the order of the Keck diffraction limit at 2.2 /xm, this 
introduces a Strehl dependence of the i^'-band flux measurements. The data show that the 
Strehl on NGC 4151 is systematically lower than on the calibrator stars, because the stars are 
significantly bluer then the red AGN and deliver more photons to the visible wavefront sensor 
(WFS; note that we matched calibrator stars and AGN in the infrared brightness to get good 
visibility accuracy). Furthermore, the WFS probably detects some extended flux from the 
host galaxy. 

To account for this differential Strehl, the raw data of the AGN are multiplied by 
1.45, before the flux calibration of the raw data against the calibrator stars. This empirical 
correction factor for the differential Strehl (between stars and AGN) of 1.45 appears to be 
stable over time, within the precision of the flux count estimation. It fits to the data of the 
two campaigns targeting NGC 4151, where qua s i-simul taneo us K-hand photorn etry from a 
separate telescope is available: iKishimoto et al.l ( l2009bl ) and lSwain et al.l (120031 ). While the 
former publish quasi-contemporaneous UKIRT photometry, we reanalysed the archived KI- 
data of the latter paper, and compared the Fatcat photometry ag a inst t he simultaneous 
Magnum single telescope photometry, as presented by lKoshida et al.l ( l2009l ). To cross-check 
that by applying the correction factor, we account properly for the loss of light due to the 
spatial filtering of our i^'-band data, we also flux-calibrated the spatially unfiltered 1.6 /im 
(if)-band flux of the tip-tilt imager next to Fatcat, which feeds a tip-tilt stabilizing loop 
during KI measurements. Both this if - band flux (50 mJy), and the corrected i^-band flux 
(95 mJy) fit the corresponding values in IKishimoto et al.l (l2009bl ). with the same flux offset 
of a factor of 1.9, which we account to intrinsic flux variability. In Fig. [21 we used for the 
i^-band photometry of the three visibility measurements the corrected Fatcat photometry 
for the Swain et al., and our data, and the UKIRT photometry from Kishimoto et al. 2009, 



http : //nexsci . caltech . edu/ sof twaxe/KISupport/v2/v2sensitivity . html 
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together with the respective photometric uncertainties. The flux contanainatioii of the host 
galaxy in the i^'-band is with about 1 % neglectable (IKishimoto et aLll2009bl ). The flux 
correction factor of 1.45 for the Fatcat i^-band data of NGC 4151 seems to be stable over 
the years, and therefore not too sensitive to the actual observing conditions, as long as a 
similar calibration scheme is used. This is probably due to the fact, that in case of poor, 
unstable AO performance the interferometric measurement would not work efficiently. But 
the correction factor should not be applied to future datasets without careful checking of the 
effective differential Strehl between calibrator and target during the respective observation. 
In particular, the correction factor should be re-estimated for KI observations of different 
targets. 



3. Discussion of the results 



The calibrated visib il ities are plotted i n Fig. [T] together with previously published data 
from ISwain et al.l (120031 ) ; IKishimoto et al.l (j2009bl ) . It is apparent that the different inter- 
ferometry campaigns measured very similar visibilities, despite of the significantly different 
total flux measured (Fig. |2]). By fitting models, we focus on three different properties of the 
size constraints which can be derived from the visibilities: the amount of unresolved fiux, 
the size of the resolved fiux emission structure, and the deprojection of this size. Our models 
consist of a variable fractional amount of (unresolved) nuclear emission [FR^^^ = — 1) 
along with resolved (overresolved) emission at (beyond) centro-nuclear radii -Rdust (Fig-E]). 
We simply refer to the extended component as dust, being the most likely origin of it. 



3.1. Modeling the dust emission 

Model 1: Already a single visibility below unity shows that some fraction of the fiux is 
resolved by the interferometer. This corresponds to FR^^^ < 1. The most compact bright- 
ness distribution, still consistent with the data, would put FR^^^ ~ 90 % of the fiux at radii 
smaller than 0.35 mas (< 30 mpc) below the resolution limit of the KI ("point source"), 
and the remainder of the fiux outside of -Rdust ~ 400 mpc to be completely resolved out. 
Since radii of about 30 mpc are significantly smaller than the dust subli mation radius in 



NGC 4151 , even for large graphite grains (see the recent discussions in IKishimoto et al. 



20071 : iNenkova et al.l l2008bl ) , this scenario would suggest that the i^'-band emission is dom- 
inated by hot, and potentially partially non-thermal nuclear emission from the innermost 
accretion disk zone. 
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Fig. 1. — Calibrated visibility data from our observations (blac k), and from two other mea- 
surements (blue, red: ISwain et al.ll2003t iKishimoto et al.ll2009bl respectively), plotted versus 
the projected baseline length. The dotted errorbars show the individual data points, with 
their estimated uncertainties. The three larger, solid errorbars show the mean V^^-result 
of each campaign, and the adequate accuracy bias of at least 0.03. The smaller statistical 
errors reflect the differential precision of the observations within a night, the larger accuracy 
bias needs to be taken into account if absolute size scales are fitted, and if observations from 
different nights are compared, minimized face-on ring models are overplotted, the best 
fit diameters are given. A compact flux contribution from the AGN on the order of up to 
25 % does not alter the results significantly (see Table [3] for various model fit results). The 
total K-hand flux shown comprises the flux contributions from the accretion disk and from 
the circum-nuclear dust (see Sect. [2] for flux calibration details). 
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NGC 4151 : size-luminosity relation 
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Fig. 2. — Dust ring radii, derived from fitting a face-on ring to the calibrated visibility data. 
Color-coding matches Fig. [T] In the lower left, a L*^'^-dependence is shown for comparison. 
Note, that the plotted radii neither take into account projection effects nor variable nuclear 
flux contributions. An inclined dust ring would have larger centronuclear radii. A nuclear 
flux contribution, varying between 10 - 25%, could explain the apparent slope of the three 
datapoints, as discussed in Section 13.21 
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Geometric constraints from the visibilities of NGC 4151 

Resolved emission (ring) + nucleus 




>400mpc 

unresolved nucleus 

« R„,„ < 30 mpc 

Fig. 3. — Sketches (not to scale) of the geometric constraints on the extended i^'-band 
emission of NGC 4151, seen through the 50 mas spatial filter of the Keck interferometer. 
Details are given in Section 13.11 Note, that the overresolved emission in model 1 might be 
diffuse, or more spherically than drawn, but the emission needs to origin from centronuclear 
radii larger than 400 mpc to be completely resolved out, and to not contribute to the V^- 
budget. 
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But the derivation of FR^^^ from a few visibilities is degenerate. Reducing FR^^^ re- 
quires only to shrink -Rdust to stay consistent with the data. The visibility alone does not 
constrain FR^^^ any further than being < 0.9. However, recent spectro-photometric and 
NIR-polarimetric data appear to limit FR^^^ to less than 25 %. In these studies, the 
central i^-band flux appears dominated by moderately warm thermal emission (~ 1500 K), 
likely from circum-nuclear dust clouds (see the detailed discussions and references of the ob- 
servation al support of dust being the dominating origin of the K-ha nd emission in Seyfert 1 
nuclei in iBarvainisI 119871 : iRiffel et al.l l2009t iKishimoto et al.ll2009al Jbl). Since our KI pho- 
tometry and visibilities from the Dec 2008 run are consistent with these earlier findings, we 
assume in the following FR^^^ < 0.25, which means that both total and correlated flux are 
dominated by extended emission with an emission size similar to the angular resolution of the 
interferometer. Such resolved emission would result in decreasing visibilities with increasing 
baselinelength. Indeed such a trend is marginally apparent in both new datasets with high 
intranight differential visibility precision (Dec 08, and May 09 in Fig. [1]). 

Model 2: aims at estimating the order of magnitude of -Rdust- To translate visibilities into 
size scales, a simple model of the extended emission is required. Gaussians and rings fit the 
data equally well, and give the similar -Rdust.circ of ~ 45 mpc (0.5 • FWHM f» Rring = 0.55 ± 
0.1 mas. Fig. [Hand Table [3]). This size estimate is largely unaffected by an FR^^^ < 25%, or 
by the radial extension of the ring as long as the ring radius is larger than its extension. The 
key difference between Gaussian and ring models is, that the Gaussian not only defines a size 
scale, but at the same time it constrains the emitting area. The fitted FWHM translates into 
a solid angle of 0.92 mas^. However, a brightness temperature calculation rejects a smooth 
Gaussian-like emission surface profile. Assuming a black body temperature of 1300 K, this 
surface would result in about 500 mJy. This is about a factor of 2.5-10 larger than the 
measured dusty i^-band emission of NGC 4151, depending on its status of nuclear activity 
(Fig. 12]). If the dust emission is optically thick, this discrepancy means that the radiating 
surface area is significantly smaller than the solid angle of the fitted Gaussian. If we were 
to shrink the FWHM of the Gaussian to match the observed fluxes, the resulting increased 
visibilities would overshoot over the data. Higher dust temperatures would strengthen this 
argument against a smooth Gaussian brightness distribution. T hus we conclude t hat the 
black body component, apparent in the spectro-photometry (e.g. iRiffel et al.l 120091 ) . is not 



produced by an approximately spherical, smooth dust distribution (which would resemble 
a Gaussian in projection). The planar geometry of the ring model seems more reasonable. 
This matches the expectation for observing a type 1 AGN. Optically thick clouds, as used 
in modern torus models, would reduce the radiating surface as well (without changing the 
-Rdust estimate). The ring morpholog y also appears re asonable due to the lack of significant 
extinction towards the very nucleus (ILacy et al.lll982[ ). 



- 12 - 



Model 3: The -Rdust,circ estimate from the previous paragraph assumes a face-on, circular 
dust distribution. However, given the single KI baseline, the size constraints of Kl-data are 
primarily one-dimensional along the position angle of the baseline (about 40°). The resolved 
source could be substantially larger along the orthogonal (SE-NW) direction. In fact, the 
KI baseline position a n gle roughly co-aligns wit h the ionization cone axis of symmetry of 
NGC 4151. iDas et al.l (120051 ): iRiffel et al.l (120091 ) report position, opening and line-of-sight 
inclination angles of the cone to be ~ 60°, ~ 70°, and ~ 45°. If the circum-nuclear 
dust, seen as extended component by the interferometer, is responsible for the formation 
of the ionization cone, the opening angle would translate into a torus half-height of 0.7 
times the inner radius. Such a geometry matches the thin ring model, if the extended K- 
band emission is dominated by the hottest dust exposed to direct nuclear illumination. 
If KI and the projected minor axis were exactly co-aligned, the size estimate from the 
previous face-on model (-Rdust,circ ~ 45 mpc) needs to be corrected for the inchnation angle to 
-Rdust.inci ~ 60 mpc (see Table [3]). The position angle of the inner radio-jet, a possibly better 
constraint on the AG N axis of symmetry, lies with a PA of ~ 80° within this ionization cone 
(IMundell et al.ll2003l ). but slightly off the KI position angle. Summarizing, the intrinsic 3d 
inner torus radius -Rdust is constrained to -Rdust,circ < -Rdust < -Rdust,mci- Taking into account 
the orientation constraints of radio jet and ionization cone, we estimate -Rdust = 55 ± 5 mpc. 
The significant inclination of the nuclear axis of symmetry, indicated by the inclination of 
the ionization cone, is supported by the tirn e- variable, spectral classification of NGC 4151 
as Seyfert 1.5 - 1.8 (IShapovalova et al.ll2008l ). 



3.2. On the luminosity dependence of the torus size 



In general, both y-to-K-continuum reverberation and NIR interferometric measure- 
ments of nearby AGN appear to confirm a L °'^-dependence of the extended -ft^-band emissio n 
size on the nuclear illumination luminosity (ISuganuma et al.ll2006l : iKishimoto et al.ll2009bl ). 
In addition, the NIR SED bump of AGN can be fitted consistently with a black body 
profile of about 1500 K , for AGN samples covering orders of magnitude in nuclear power 
jKobavashi et al.lll993h . These are strong arguments for the i^-band being dominated by 
thermal dust emission close to the sublimation limit. The strong nuclear emission variabil- 
ity of NGC 4151 by several 100 % on yearly timescales makes it an ideal target to do the 
next step, and to study the relation between the nuclear emission and the dust emission 
size in detail. Sequential dust reverberation studies in several AGN did find changes of the 
response delay time, suggesting a cha nge of -Rdust in reaction to different nuclear illumination 
(lBarvainislll992l : lKoshida et al.ll2009l ). However, given a -Rdust of several tens of mpc, it would 
require unreasonable fast dust bulk motion (at 0.1 c) to follow the luminosity changes ac- 
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cording to a L'^'^-law. Therefore, cyclic dust grain des truction and for mation was suggested 
to explain both the observed infrared flux variations (IBarvainisI Il992l ) and the large scatter 
and relatively poor L ^'^-dependence in the dust reverberation size estimates of NGC 4151 
jKoshida et aliboool ). 



For the first time, we can compare these claims of a changing -Rdust against a multi-epoch 
visibility dataset, taken at different brightness states (Fig. [2]) of NGC 4151. Assuming a 
simple LP-^ scaling, the flux variation between the two new datasets (-Fk.Mayog / -pK.Decos = 
1.9 ± 0.3) would result in a size increase by a factor of 1.4 ± 0.1. Such large size variations do 
not fit the multi-epoch dataset, since the typical precision of the centronuclear radius estimate 
is < 10 %. In addition, the flux- dependence of -Rdust is constrained to i? oc ]_,{-'^-'^^-'^-^) by 
the data. Thus, if indeed we would observe a variation of -Rdust, then the trend would predict 
larger sizes at lower luminosity states, in contrast to a -L°'^-dependence. 

However, the apparent slope can be explained by -F-R^c variability alone, witho ut the 
need for a variation of -Rdust between the three observations. iKishimoto et al.l ( 120071 ) mea- 
sured a significantly higher nuclear flux contribution to the -fT-band at a brighter state of the 
nucleus than during the observation of iRiffel et al.l (120091 ) (20-25% with respect to 10%), 
although both teams fitted in comparable ways very similar black body temperatures to 
the NIR-emission. The comparison of the two new precise visibility measurements reveal 
a small ^^-offset of 0.04 at the same angular resolution (compare the black and red data 
in Fig. [TJ Such a difference would result from changing FR^^^ between both measurements 
from a 10 % level in Dec 08 to the 20 % level in May 09. The small difference, slightly 
larger than the conservative internight accuracy estimate of 0.03, can be accommodated by 
a change of the compact FR^^^, without further assumptions on the structure of the ex- 
tended emission, in particular without changing the morphology of the extended brightness 
distribution. However, a change of the thermal flux by a factor of two, without changing the 
size of the radiating surface, would require a temperature change by ~ 1.2. As far as we 
know, such a luminosity dependent dust temperature variation has not been observed yet in 
NGC 4151. In contrast, an increased radiation surface of a cloudy torus of fixed size would 
point to larger cloud covering factors or an increased number of efficiently heated clouds in 
brighter states. 

Optical polarization, and NIR-continuum variations appear correla ted with visible nu- 
clear continuum variations and enable reverberation experiments (see iGaskell et al.l l2007l . 
for a polarization reverberation experiment). Therefore, a variation of the fractional FR^^^ 
cannot be explained by a varying line-of-sight extinction towards the nucleus (e.g. due to 
passing clouds), because such an apparent nuclear luminosity change should not result in 
the observed correlated change of optical polarization, and NIR-continuum fluxes. However, 
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in a phase of nuclear brightening, the torus response might lag the nuclear flux rise by the 
reverberation delay. This could explain a momentary increase in FR^^^ and in V^, even if 
the intrinsic -Rdust does not change. 

50mpc equa ls 60 days of li g ht tra vel, which matches the first four dust reverberation 
measurements of iKoshida et al.l (120091 ). Because these earlier dust reverberation measure- 
ments span already a factor of six in optical luminosity, they also su ggest a luminosity- 
variation independent dust emission size. The first KI measurement by ISwain et al.l (j2003[ ) 
is quasi-contemporaneous with one of these dust reverberation measurements. Thus both 
methods appear to agree very well. However, the former authors report on four subsequent, 
significantly smaller reverberation size estimates, the shortest being ~ 30 mpc only. Rever- 
berat ion size estimates are biased to the innermost responding material ( IGaskell fc Sparke 
19861 ). Furthermore, if the responding structure is rather flat and inclined with respect to 
the line-of- sight, as expected for NGC 4151, the first response can set in considerable faster 
than light travel time corresponding to -Rdust- However, in case of an isotropic illumination, 
and a fixed dust geometry, such geometric effects could only explain a fix offset between 
reverberation and interferometric radii. 

This discussion implies that dust reverberation sizes may often underestimate the loca- 
tion of the bulk of the circum- nuclear dust, which dominates the i^'-band emission. Indeed, 
if we use our -Rdust for NGC 4151, at lea s t three of the four interferometric -Rdust estimates of 
the AGN, observed by lKishimoto et al.l ( l2009bl ). exceed the L*^'^-fit of the reverberation size 
by factors of about 2 ± 0.5. The four th object, Mrk 23 1, likely has an inclined geometry, 
and larger, de-projected -Rdust as well ( ISmith et al.l 120041 ). The direct interferometric -Rdust 
estimate appears as the more robust measurement of the average location of the bulk of the 
warm circum-nuclear dust (if projection ambiguities are well constrained). In contrast, the 
smallest reverberation sizes are smaller than the interferometric sizes, and the reverberation 
sizes vary beyond the scatter of the interferometric size estimates. We suggest, that the 
reverberation sizes are probably tracing the location of momentary illuminated innermost, 
very hot dust clouds, instead of the bulk of the -ft^-band dust. If ci rcum-nuclear dust origi- 
nates in outflows from the central engine (lElitzur fc Shlosmanll2006l . and references therein), 
then the shortest dust reverberation times could derive from optically thick dust in such 
outflows inside -Rdust- 



4. Conclusions 



The results of a new interferometric NIR observing campaign of NGC 4151 are pre- 
sented. Although this time in a 2 times fainter state, we managed to re-observe the AGN 
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with the Keck Interferometer thanks to recent sensitivity improvements of the instrument. 
The interferometric visibihties and size estimates were compared to previously pubhshed 
intcrfcrometric and single telescope data. The current KI datasets suggest that the major 
part of the radiation originates in a (possibly inclined) toroidal structure of an intrinsic, 
deprojected radius of about 55 ± 5 mpc, comparable to \/-to--ft'-continuum reverberation 
measurements, without constraining the morphology in further detail. 

Our datasct and its comparison to published data enables for the first time to study the 
sensitivity and response of near-infrared intcrfcrometric visibilities of an AGN to intrinsic fiux 
variations. The observations show that NlR-interferometry on bright AGN {K ~ 10 mag) 
are now feasible at a high precision of a few percent. We did not detect a significant visibility 
dependence on doubling the X-band emission of NGC 4151. This supports the notion, that, 
for an individual AGN, the size of its circum-nuclear NIR emission structure does not strictly 
depend on the respective momentary nuclear luminosity. The direct interferometric size 
estimates appear to be more robust estimates of the average location of the dust than NIR 
continuum reverberation estimates. Being fiux-independent, this average location of the hot 
dust probably relates to the sublimation radii of the AGN at its high activity state. Future 
dust reverberation campaigns, contemporaneous to direct interferometric measurements are 
highly desirable to study the circum-nuclear NIR emission region, to correctly interpret the 
reverberation measurements, and to investigate apparent changes in the dust illumination, 
heating efficiency, or covering factors. 

We are grateful to the excellent KI team at WMKO and NExScI for making these 
observations a success. K. Meisenheimer and L. Burtscher contributed helpful discussions. 
The data presented herein were obtained at the W.M. Keck Observatory, which is operated 
as a scientific partnership among the California Institute of Technology, the University of 
California and the National Aeronautics and Space Administration. The Observatory was 
made possible by the generous financial support of the W.M. Keck Foundation. The authors 
wish to recognize and acknowledge the very significant cultural role and reverence that the 
summit of Mauna Kea has always had within the indigenous Hawaiian community. We are 
most fortunate to have the opportunity to conduct observations from this mountain. The 
Keck Interferometer is funded by the National Aeronautics and Space Administration as 
part of its Navigator program. This work has made use of services produced by the NASA 
Exoplanet Science Institute at the Cahfornia Institute of Technology. This research has 
made use of the SIMBAD database, operated at CDS, Strasbourg, Prance. This research 
has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract with the National 
Aeronautics and Space Administration. 
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Table 3: Model fit parameters consistent with the data. The varying fit results of the ring 
models #2 and #3 reflect projection effects and represent general lower and upper limits on 
the intrinsic -Rdust (Sect. |3TT]) . 



Model FR^^^ -Rdust Type of dust distribution 







[mas] 


[mpc] ^""^ 




#1 


0.9 


> 5 


> 400 


overresolved 


#2 


0.25 


0.55 ± 0.05 


45 


resolved ring 


#2 





0.6 ± 0.05 


50 


resolved ring 


#2 


0.25 


0.45 ± 0.05 


35 


resolved Gaussian 


#2 





0.5 ± 0.05 


40 


resolved Gaussian 


#3 





0.75 ± 0.05 


60 


resolved ring, 45° inclined, 
minor axis co-aligned with the KI baseline 



"from Ned: Cosmology-corrected redshift of 0.00414 against the 3K microwave background refer- 
ence frame results in angular-size distance of 16.9 Mpc, and 82 mpc/mas, respectively, with Hq = 
73kmsec^^ Mpc^^, riniattcr — 0.27 



